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Abstract: 

We demonstrate a kind of linear superposition for a large number of nonlinear equations, both contin- 
uum and discrete. In particular, we show that whenever a nonlinear equation admits solutions in terms 
of Jacobi elliptic functions cn(x,m) and dn(x,m), then it also admits solutions in terms of their sum as 
well as difference, i.e. dn(x,m) ± ^/m cn(x , m) . Further, we also show that whenever a nonlinear equation 
admits a solution in terms of dn 2 (x, m), it also has solutions in terms of dn 2 (x, m) ±^/mcn(x, m) dn(x, m) 
even though cn(x, m) dn(x, m) is not a solution of that nonlinear equation. Finally, we obtain similar 
superposed solutions in coupled theories. 
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Introduction: Linear superposition principle is one of the hallmarks of linear theories which does 
not hold good in nonlinear theories because of the nonlinear term. For example, even if two solutions are 
known for a nonlinear theory, their superposition is in general not a solution of the nonlinear theory. The 
purpose of this letter is to point out a kind of superposition for a large number of nonlinear equations. 
In particular, there are several nonlinear equations, both discrete and continuum, which are known to 
admit exact periodic solutions in terms of Jacobi elliptic functions cn(z, m) as well as dn(x,m), where 
m denotes the modulus of the elliptic function [T]. Many of these solutions have found application in 
several areas of physics [2J|3]. In particular, we have examined a large number of nonlinear equations, both 
continuum and discrete, which admit both cn(x, m) and dn(x, to) solutions and find that in all these cases 
even dn(x,m) ± y / mcn(x,m) is also an exact solution. We have also examined a number of coupled field 
theories and have obtained superposed coupled solutions of the form dn ± y^mcn in both the fields. 

Further, we have also examined a number of continuum field theories which admit dn 2 (x,m) as a 
solution and find that such theories also admit dn 2 (x, m) ± ^ [ /mcn(x, m) dn(x, to) as solutions. While this 
cannot be treated as a linear superposition of two solutions (since cn(x, m)dn(x, to) is not a solution of 
these models), we find it rather remarkable that such solutions exist, without exception, in a number of 
continuum field theories (including several coupled models) that we have looked at so far. 

In this letter we only discuss a few selected examples, several other examples will be given elsewhere 
in a longer version. To begin with, we discuss one continuum (quadratic-cubic nonlinear Schrodinger 
equation or QCNLS) and one discrete model (saturated discrete nonlinear Schrodinger equation or DNLS) 
both of which admit cn(x,m) as well as dn(x,m) as solutions and show that these models also admit 
dn(x,m) ± y/m cn(x, m) as solutions. Note that stable kink [2] and chaotic soliton solutions under nonlin- 
earity management [5] are known for the QCNLS equation which arises in such diverse physical systems 
as nonlinear optics, chemical kinetics, matter-radiation interactions and mathematical ecology [HE]. The 
DNLS equation arises in the context of optical fibre communication [6]. We then discuss a coupled NLS- 
MKdV model and show that it has dn(x, to) ± ^/m cn(x, m) type coupled solutions in both the fields. Here 
MKdV refers to the modified KdV equation [2]. 

Subsequently, we discuss the Korteweg-de Vries (KdV) equation which is known to admit dn 2 (x,m) as 
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a solution [2j [3] and show that this model also admits dn 2 (x, m) ± ^/mcn(x, m) dn(x, m) as solutions even 
though cn(x, m) dn(x, m) is not a solution of the KdV equation. We also show that the coupled NLS-MKdV 
model, not only admits dn(x,m) ± T/rncn(x,m) but also dn 2 (x,m) ± - v /mcn(x, m)dn(x, m) type solutions 
in both the fields. Further, we show that a coupled NLS-KdV model has superposed solutions of the form 
dn(x,m) ± y / mcn(x,m) - dn 2 (x,m) ± \/mcn(x,m) dn(i,m). Finally, we briefly mention possible reasons 
why such a linear superposition works in several nonlinear theories, 
dn ± y/rncn as Exact Solutions 

We first discuss the quadratic-cubic nonlinear Schrodinger equation (QCNLSE) and then the discrete 
NLS equation (DNLSE) and show the existence of such superposed solutions in both the cases. 
QCNLSE: We start with the quadratic-cubic nonlinear Schrodinger equation 

iut + u xx +gx\u\u + g 2 \u\ 2 u = 0. (1) 

One of the exact moving periodic solutions to this equation is 

u = (Adn[f3(x -vt + <5i), m] + B) exp[-i{ut - kx + 5)} , (2) 

provided 

g 2 A 2 = 2/3 2 , g 2 B 2 = (2-m)/3 2 , gi = -3Bg 2 , u = k 2 + 2(2 - m)f3 2 , v = 2k . (3) 

Here 8, S\ are two arbitrary constants. 

Similarly, another exact moving periodic solution to the QCNLSE, Eq. (TjQ), is 

u = {Ay/mcn[f3{x — vt + 5\), m] + B) exp[— i{uit — kx + 5)], (4) 

provided 

g 2 A 2 = 2(3 2 , g 2 B 2 = (2m- l)/3 2 , 9l = -3Bg 2 uj = k 2 + 2(2m - 1)(3 2 , v = 2k . (5) 

Notice that the cn(x,m) solution only exists if 1/2 < m < 1. We now show that, remarkably, even though 
cn(x, m) solution does not exist if m < 1/2, a linear superposition of cn(x, m) and dn(x, m) is still an exact 
solution over the entire range < m < 1, i.e. 

u = ^— dn[/3(x — vt + <5i), m] + — \fm cn[/3(x — vt + 5i), m\ + B^j ex.p[—i(ujt — kx + 5 2 )] , (6) 
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is an exact solution to the QCNLSE, Eq. ([TJ, provided 

D = ±A, g 2 A 2 = 2(3 2 , g 2 B 2 = (l + m)/2p 2 , 

9l = -3Bg 2 , oj = k 2 + (l + m)/3 2 , v = 2k. (7) 

Here cn(x,m) and sn(x,m) are periodic functions with period 4K(m), dn(x,m) is a periodic function 
with period 2K(m), with K(m) being the complete elliptic integral of the first kind [T]. It is worth noting 
that the frequency oj of the three solutions (i.e. cn, dn and dn ± y'm cn) is different except at m = 1. We 
thus have two new periodic solutions of QCNLSE depending on if D = A or D = —A. 

Few remarks are in order here which are in fact valid for all the models (both continuum and discrete) 
that admit such solutions. 

1. Both the solutions dn + y^cn and dn — ^/rn cn exist for the same values of the parameters. 

2. In the limit m = 1, the three solutions dn, cn as well as dn + ^/mcn go over to the well known pulse 
(i.e. sech) solution. 

3. In all the continuum models admitting such solutions, the factors of 2 — m or 2m — 1 which appear 
in the dn and cn solutions, get replaced by the factor of (1 + m)/2 in the dn ± ^fm cn solutions. 

4. On the other hand, in the discrete theories admitting such solutions (see below), the factor of dn(/3, m) 
or cn(/3,m) appearing in dn and cn solutions, gets replaced by a factor of [dn(/3,m) +cn(/3,m)]/2 in 
the dn ± -^/mcn solutions. 

5. In view of the above two points, either the frequency oj or the velocity v (or could be even both) are 
different for dn, cn and dn ± -y/mcn solutions except at m = 1. 

Saturated Discrete NLS Equation: We now consider a discrete nonlinear equation which admits 
both dn and cn solutions and show that it also has superposed solutions dn ± ^/mcn. In particular, 
we consider saturated DNLS equation which has received great attention in the context of optical fibre 
communication [6] 

/ i r "i I *^^"t~h t / \ 

idUn/dt + [U n+1 + U n -l\ + — — : prU n = . (8) 
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One well known periodic solution to this equation is [7] 

u n = Ada\J3(n + S t ), m\e^ t+ ^ , (9) 

provided 

■ 9 9,„ s „ 2K(m) dnf B,m) , . 

A 2 cs 2 /3,m =1, p = —±-L, u = -u = -2 ^' ; 10 

iVp cn^ (p,m) 

The other solution is 

u n = AV^cn[/3(n + 5i), m]e~ i(a; * +<5) , (11) 



provided 

.9, . „ 4K(m) cn(8,m) 

A 2 d S 2 (f3,m) = l, (3 = — ±-L, u = -v = -2 £ \ - ( 12 ) 

iV p dn (p,m) 



Remarkably, even a linear superposition of the two is also an exact solution to the saturable DNLS Eq. 
([8]), i.e. it is easily shown using the recently derived identities for the Jacobi elliptic functions [8] that 

u n = (-dn[/3(n + 6i),m] + jV^cn[P(n + S 1 ),m]\- i(iJt+S \ (13) 

is also an exact solution to Eq. © provided 

A-K(m\ 4 
B = ±A, ra ) + ds(ft m)f = 4 , = , „ = -„ = - dn(ft ,„ )+cn(ftm) ■ (14) 

Here iV p is the spatial period of the system, cs(/3, m)=cn(/3, m)/sn(/3, m) and ds(/3, m)=dn(/3, m)/sn(/3, m). 
As remarked earlier, observe that if we replace dn(/3, m) and cn(/3, m) by [dn(/3, m) + cn/3, m)]/2 in relations 
(|10p and (|12p . we recover relations (|14p . Further, the frequency u; of the three solutions dn, cn and 
dn db y/mcn is different except at m = 1. 

Coupled NLS-MKdV Model: We now show that the coupled NLS-MKdV model admits dn^t^/mcn 
type solution in both the fields. 

In this case the field equations are given by 

i 1 9 9 

iut + u xx + g\u\ u + auv =0, 

v t + v xxx + 6v 2 v x + <yv(\u\ 2 ) x = . (15) 



5 



These coupled equations admit four solutions with u being either en or dn (multiplied by an exponential) 
and similarly v can be either cn or dn. For example, one of the solutions is 

u(x, t) = Aexjp[—i(u>t — kx + <5)]dn[/3(x — ct + 5i),m] , 

v(x, t) = By/mcn\]3(x — ct + 6i),m] , (16) 

provided 

c = 2k = (2m- l)/3 2 , a = k 2 - (2 - m)/3 2 + (1 - m)aB 2 , (17) 
gA 2 + aB 2 = 2(3 2 , -/A 2 + 3B 2 = 3(3 2 . (18) 
Remarkably, even a linear superposition 

u(x, t) = — exp[— i(ut — kx + 5)] ( Adn[f3(x — ct + Si), m] 



+Dy/mcn[j3{x — ct + Si),m]j , 

v(x,t) = ^(bAd\]3{x - ct + 5i),m] + F^/mcn[f3(x - ct + 5{),m]j , (19) 
is an exact solution of coupled Eqs. (|15p provided 

c = 2k = (1 + m)(3 2 /2 , lo = k(k - 2) , D = ±A, F = ±B , (20) 



while A, B are given by Eq. (|18|) . Note that the signs of D = ±A and F = ±i? are correlated. 
dn 2 =t ^/mcn dn superposed Solutions 

We now discuss one example of a continuum field theory which admits dn 2 as a solution and show that 
the same model also admits dn 2 ± y'mcndn as solutions, even though cndn is not a solution of the model. 
KdV Equation: It is well known that one of the exact solutions to the KdV equation [2j 

ut + Uxxx + guux = , (21) 

is 

u = Adn 2 [f3(x-vt + 5i),m], (22) 

provided 

gA = l2/3 2 , v = 4(2- m)/3 2 . (23) 



Remarkably, even a linear superposition, i.e. 

u = i^4dn 2 [/3(x -vt + 8i),m] + B^cn[f3(x - vt + 61), m] dn[/3(x -vt + 6i),m]j , (24) 

is an exact solution to the KdV Eq. (|2"T|) provided 

B = ±A, gA = 12/3 2 , v = (5-m)/3 2 . (25) 

We thus have two new periodic solutions of KdV Eq. (|2"Tj) depending on if B = A or B = —A. 

Several remarks are in order here which are in fact valid for all the solutions of the form dn 2 ± ^Jmcli dn. 



1. We find that all the models which admit dn 2 as a solution, also admit dn 2 ± -^/mcndn as a solution 
even though -^/mcn dn is not a solution of these models and that both of the new solutions exist for 
the same values of the parameters. 

2. In the limit m = 1, the two solutions dn 2 and dn 2 + -y/mcndn go over to the well known pulse (i.e. 
sech 2 ) solution while dn 2 — -y/mcndn solution goes over to the vacuum solution u = 0. 

3. The factors of 2 — m and y/l — m + m 2 which appear in the dn 2 solution, get replaced by the factor 
of (5 — m)/4 and vT+ 14m + m 2 /4, respectively, in the dn 2 ± \/mcndn solutions (see below). As 
a result the velocity v or frequency u (or even both) for the solutions dn 2 and dn 2 ± -^/mcndn are 
different (except at m = 1). 



We now show that the same coupled NLS-MKdV model as given by Eq. (|15p . also admits dn 2 as well 
as dn 2 ± -y/mcndn type solutions in both the fields. 

It is easily shown that the coupled system as given by Eq. (|15j) admits an exact solution 

u = ^4dn 2 [/3(x - ct + 5i),m] + F^j exp[-i(ut - kx + 5)} , 

v = Bdn 2 [p(x-ct + 6 1 ),m} + D, (26) 



provided 



a = 3/2, -fA 2 = -3B 2 , 7 = 2 5 <0, (z - y)B 2 = 2/3 2 , 

c = 2k = 4[(2 - m) + 3z]/3 2 , w = A; 2 - [4(2 - m) + 9y + 3z]/3 2 , 



I? F [-(2-m)± v / T 

5' y = I = 3 • (27) 
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Note that Adn 2 + B is not a solution of either the NLS or MKdV uncoupled field equations even though 
the coupled system admits such a solution. 
Remarkably, even a linear superposition, i.e. 

u= ( F + —dn 2 [(3 (x - ct + 5i),m] 

Q \ 

+— v / mcn[/3(x — ct + 6\), m]dn[(3(x — vt + 5i),m] J exp[— i(uot — kx + 5)} , 
v = D + — dn 2 [/3(x - ct + <5i),m] 

+— v / mcn[/3(x - ct + m]dn[/3(x - ct + ^),m] , (28) 
is an exact solution to Eq. (|15j) provided 

G = ±A, H = ±B, c = 2k = [5 - m + I2z]p 2 , 

, 9 r, \ n ^ 9 F f-(5 - m) ± + 14m + m 2 l . . 

w = A; 2 - [(5 - m) + 9y + 3z]/3 2 , y = - = ^ ; V ^ (29) 

while rest of the relations are exactly the same as those given by Eq. (|27p . Note that the signs of G = ±A 
and H = ±B are correlated. It is worth reminding once again that neither Adn 2 + B nor A^/mcndn is an 
exact solution of either NLS or MKdV field equations. 

Coupled NLS-KdV Fields: Finally, let us consider the following coupled KdV-NLS field equations 

i 1 9 

iut + u xx + g\u\ u + auv = 0, 

Vt + v xxx + 6vv x +-fv(\u\ 2 ) x = 0, (30) 

where u and v are the NLS and KdV fields, respectively. These equations admit coupled solutions of the 
form dn - dn 2 and cn - dn 2 . For example, it is easily checked that 

u(x, t) = Aexp[— i(uit — kx + <5i)]dn[/3(x — ct + 5), m] , 

v(x,t) = Bdn 2 [l3{x-ct + 5i),m] + D, (31) 



is an exact solution to the coupled Eq. (|30j) provided 

gA 2 + aB = 2f3 2 , jA 2 + 6B = 12/3 2 , (32) 
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c = 2k = 4[2 - m + 3z]/3 2 , u = k 2 - (2 - m)/3 2 - azB , z = ^ . (33) 

Remarkably, the same model also admits interesting superposed solutions of the form dn db i/mcn — 
dn 2 ± -y/mcndn. In particular, it is easily checked that 

u(x, t) = — exp[— i(uit — kx + 5)] ^Adn[j3(x — ct + 6i),m] 

+Hy/mcn[f3(x — ct + <5i), m] 

v(x, t) = — dn 2 [/3(x — ct + 5i),m] 
F 

+—^/rncx).[(3{x - ct + St), m]dn[/3(x - ct + 5i),m] + D , (34) 



is an exact solution to the coupled field equations (|30|) provided 

z = —, H = ±A, F = ±B, c = 2k = [(5 - m) + 12z]/3 2 , 
B 

w = A; 2 - ^ ^ -_[(!_ m ) + 4*]fl , (35) 

and further A, B are again given by Eq. (132p . 

Conclusion. In this letter we have shown that a kind of linear superposition holds good in the case of 
quadratic-cubic NLSE and saturated discrete NLS equations in the sense that these models not only admit 
dn(x,m), cn(x,m) but even dn(x,m) ± y / mcn(x,m) type superposed solutions. In fact, by now we have 
examined a large number of both discrete and continuum nonlinear equations such as NLS, MKdV, mixed 
KdV-MKdV system, A <fi field theory, Ablowitz-Ladik equation, saturable discrete NLSE, discrete A <j) field 
theories, among others, and in all these cases we have obtained such superposed solutions. Furthermore, 
even in several coupled theories such as NLS-MKdV we have obtained dn(x, m)±y / mcn(a;, m) type coupled 
solution in both the fields. 

Further, we have shown that for KdV equation, which is known to admit a periodic solution in terms of 
dn 2 (x, m), also admits a solution of the form dn 2 (x, m) ±y / mcn(x, m)dn(x, m). By now, we have examined 
a number of nonlinear equations like quadratic NLSE and in all these cases we have obtained such solutions. 
Further, in several coupled theories like NLS-MKdV, quadratic NLS-KdV, and others too we have obtained 
dn 2 (x,m) db y / mcn(x, m)dn(x, m) type solutions. We have also looked at the NLS-KdV coupled system 
and obtained mixed solutions of the form dn(x,m) dz y/mc^x^m) - dn 2 (x,m) dz ypm cn(x, m) dn(x, m). 



In addition, we have examined several other coupled equations such as MKdV-KdV, NLS-quadratic NLS, 
MKdV-quadratic NLS, among others, and there too we have obtained such superposed solutions. 

Many of the nonlinear equations that we have looked at have found wide application in several inter- 
esting physical situations [TJ [21 El El E] including fibre optics, chemical kinetics and mathematical ecology. 
It would be worthwhile to enquire if the new solutions that we have obtained have some physical relevance, 
e.g. in optical fibre communication and related contexts [BJ. As a first step in that direction, it is important 
to examine the (linear and nonlinear) stability of these newly found solutions in the various models. 

We would like to reemphasize that what we have obtained is only a kind of linear superposition and 
not the full superposition as obtained in the linear theories. However, we find it remarkable that in spite 
of the nonlinear terms, even a kind of linear superposition holds good in these nonlinear models. 

What could be the possible reason why such a linear superposition is possible? We surmise mathemat- 
ically that the main reason is that cn(x, m) and dn(x, m) functions are quite similar and both of them as 
well as their derivatives are identical at m = 1. This is in contrast to the sn(x, m) elliptic function which is 
different from both cn(x, m) and dn(x, m) functions at any value of m and that is why cn(x, m) ± sn(:c, m) 
superposition does not seem to work. We have checked it in the various examples mentioned above. Per- 
haps there is a deeper physical reason to all these findings which needs to be explored. 

This work was supported in part by the U.S. Department of Energy. 
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